Environmental Changes, Co-extinction, and Patterns in the Fossil Record 
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We introduce a new model for large scale evolution and extinction in which species are organized 
into food chains. The system evolves by two processes: origination/speciation and extinction. In the 
model, extinction of a given species can be due to an externally induced change in the environment 
or due to the extinction of all preys of that species (co-extinction). The model is able to reproduce 
the empirical observations, such as the statistical fractality of the fossil record or the scale-free 
distribution of extinction events, without invoking extinctions due to competition between species. 



The identification of the mechanisms responsible for 
large-scale evolution and extinction is a topic of heated 
debate [0-^O| ■ The basic problem can be summarized by 
two questions. The first one centers on the cause of mass 
extinction: Is it external to the system — e.g., due to ex- 
traterrestrial impacts — or is it internal and due 
to the nonlinear dynamics of the ecosystem |h^~|Tg|]? The 
second question centers on the mechanisms for selecting 
the species that become extinct: In standard extrapola- 
tion of Darwinian evolution theory, it is assumed com- 
petition leads to the extinction of less fit species |L5| [20j , 
but some authors argue that competition might be not 
a determinant factor for macroevolution QJ^H ■ The fos- 
sil record has yet to answer these questions unequivocally 
Til- 24 1 . In qualitative modelling, all of the mechanisms 



discussed above have been considered as the explanation 
of the patterns in the fossil record. On the other hand, 
quantitative modelling approaches [fl5-18 have consis- 
tently included competition among species as a funda- 
mental mechanism. Here, we show that a quantitative 
model that does not include competition among species 
may reproduce the empirical observations, particularly 
the statistical fractality of the fossil record P,£5|, and 
the scale-free distribution of extinction sizes ppufT 

The literature on large scale species extinction reports 
on two key empirical results. First, the probability den- 
sity that a number s of species becomes extinct during a 
given time interval decays as a power law, P(s) ~ s~ T , 
with an exponent r w 2 p7| , ^8p^ |. Second, the power 
spectrum S(f) of the time series of extinction sizes also 
decays as a power law, S(f) ~ J - ' 3 , with /3 « 1 ||, which 
implies that the sequence of extinction is long-range cor- 
related. These results impose severe constraints on the 
models attempting to describe the extinction/evolution 
process. A power law decay of the probability of extinc- 
tion sizes implies that there is no characteristic size for 
extinction events, i.e. the dynamics are scale-free and 
incidents of mass extinction are likely due to the same 
mechanisms as smaller extinction events. The hypothe- 
sis that the mass extinctions are generated by the same 
dynamics as smaller extinction events is consistent with 
the self-similarity of the fossil record JoJ . 

Quantitative models have been proposed to explain 
the patterns in the fossil record. Many are based on 



the assumption that extinction events are a consequence 
of the competition between species, i. e. the least fit 
species become extinct and are replaced by new species 
[ |T5| p(i| . These changes affect the fitness of other species 
leading to bursts of extinction of all sizes. Several of 
the models @Q self-organize into a critical state in 
which many quantities are known to scale as a power law 
||18|] . However, recently it has been shown that mecha- 
nisms other than self-organized criticality, such as coher- 
ent noise can lead to power law scaling without 
requiring the system to be in a critical state. 

In this Letter, we test the hypothesis that competi- 
tion between species is not a fundamental ingredient for 
the explanation of the fossil record. This hypothesis is 
in agreement with statements that Darwinian competi- 
tion while important at the level of individuals within a 
population (microevolution) might not be relevant at the 
level of stable species (macroevolution) Thus, we 

propose a quantitative model for large scale extinction 
and evolution that does not include competition between 
species but assumes instead that the relevant mechanisms 
for macroevolution are (random) changes in the environ- 
ment |l7J,[]7J], and co-extinctions [^6| due to the interac- 
tions between species along food chains |l^,^,^8|. The 
model is able to reproduce both the power law distri- 
bution of extinction sizes and the fractality of the fossil 
record [Figs. These results suggest that competition 
between species might not be a fundamental ingredient 
for the description of the fossil record. 

The model is defined as follows. Species can occupy 
niches in a model ecosystem with L trophic levels in the 
food chain, and N niches in each level. Species from the 
first level, i = are assumed to be autotrophic (i.e., they 
produce their food through, e.g., photosynthesis), while 
species from levels I > are assumed to be heterotrophic. 
That is, a species occupying a niche in level I > feeds 
from at most k species occupying niches in level t — 1 
[Fig. [ij. We do not consider in the model any kind of 
structure of the niches within a given trophic level, that 
is, niches i and i + 1 in level I do not need to be occupied 
by similar species or to be geographically close. Finally, 
we assume that the preys of a new species are chosen 
at random from existing species in the trophic level be- 
low. The model starts with N species in level I = and 
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evolves according to the following rules: 

• Origination: every existing species gives rise, at a 
rate /i, to the creation of a new "potential" species 
that tries to occupy a randomly selected niche in 
the same trophic level or in one of the two neigh- 
boring levels. This speciation event occurs if the 
selected niche is not yet occupied by an existing 
species. Preys for the new species are selected at 
random from existing species in the trophic level 
below. 

• Extinction: at rate 1 (in some arbitrary time unit), 
a fraction p of species in the first level are ran- 
domly selected for extinction. Then, any species in 
the second level for which all preys became extinct 
also becomes extinct. This procedure is repeated 
up to level L. 

These rules imply that the number of species in the sys- 
tem is not kept constant. In particular, if the origination 
rate is smaller than a threshold value, then all species 
become extinct, i.e., the model has absorbing states [ p9[ . 
The rules for speciation imply that the origination rate of 
new successful species is proportional to the number N s 
of species in the system (leading to exponential growth, 
in agreement with the results of (|]), and to the number 
of empty niches NL — N s (which takes into consideration 
the limited resources of the system J2(J). Although the 
finite size of the system introduces competition for the 
creation of new species, the model does not involve any 
competition between existing species. 

Figure ^ shows our results for the time sequence of ex- 
tinction and origination events. The first interesting ob- 
servation is that both signals are intermittent with very 
large events appearing at a high rate. Furthermore, there 
is a strong correlation between the extinction and origi- 
nation curves, which is in qualitative agreement with em- 
pirical observations |§|3(]]. Finally, we find that the size 
of the extinction events has a distribution which decays 
with a power law tail with an exponent r = 1.97 ± 0.05, 
in agreement with empirical observations |p^7| , |l8| , ^3| . 

Next, we study the fluctuations in the number of 
species in the system [Fig. || . We find that these fluctu- 
ations are self-affinc as demonstrated by its power 
spectrum which scales as a power law. This result is 
in agreement with the perceived fractality of the fossil 
record P,§5|- 

In order to demonstrate the ability of our model to re- 
produce quantitatively the empirical data on extinction 
and origination, we compare in detail our results with the 
recent results of Ref. || . We therefore study the temporal 
correlations of extinction events for the model and com- 
pare our results with the analysis of the fossil record . 
Figure |] shows that the model results agree well with the 
empirical data, when we consider model sequences of the 
same lengths as available in the fossil record. This agree- 
ment is found for the power spectrum as well as for the 
method of detrended fluctuation analysis, which allows 



accurate estimates of correlation exponents independent 
of local trends |Q. Note, however, that once we con- 
sider longer records generated by the model, we find that 
the results crossover to uncorrelated behavior. In fact, 
the analysis of local slopes [see inset of Fig. indicates 
a similar trend for the empirical data as well, suggesting 
that extinction events might become uncorrelated at long 
time scales. 

The model proposed here is able to reproduce key sta- 
tistical properties of the fossil record, both for the ex- 
tinction and the origination of species. In contrast with 
many models in the literature, these results are obtained 
without having to assume that species have an intrinsic 
fitness, and that less fit species become extinct due to 
competition between species. In the model, mass extinc- 
tions are due to the amplification effect of predator-prey 
interactions that propagate along the food chain |l(J . In 
this framework, the extinction of some key species (due 
to environmental changes) can lead to catastrophic ex- 
tinction events. 
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the DGF for financial support. 
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FIG. 1. Schematic definition of the model. The evolution 
of the system takes place in a lattice in which each site rep- 
resents a niche in the "ecosystem". The system is organized 
into "trophic levels", a species in level £ feeds from at most 
k species in level £ — 1, except for species at the first level 
which are autotrophic. In most of the simulations there are 6 
levels with 1000 niches per level. The state of the system is 
fully described by stating the niches which are occupied by a 
specie with the list of its preys. We start the simulations with 
N ~ 50 species occupying niches in the first trophic level of 
the food chain. In the figure, the dark cells are occupied by 
a specie; the lines emerging from a cell link the species to 
its preys. The system evolves through two processes, orig- 
ination and extinction. Origination: A niche in level £ is 
randomly selected, and if a species exist there, a speciation 
is attempted: A new niche is then randomly selected in one 
of the levels £ — 1, £, or £ + 1, and if no species occupies that 
niche, a new species is created. Extinction: A fraction p of 
species in the first level are randomly selected for extinction. 
Then we remove for all species in the second level links to 
preys in the first level that have become extinct. Whenever 
all links have been removed for a species in the second level, 
it becomes extinct as well. This procedure is repeated up the 
food chain until the top level is reached. If, for the config- 
uration in the figure, the leftmost species in the lowest level 
would become extinct, then the leftmost species in the other 
levels would also become extinct. 
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FIG. 2. Extinction events are scale-free, a Time sequence 
of extinction events for the model. The lower line shows the 
individual events, while the upper curve shows the number of 
extinctions over a period of 512 time steps. Note that events 
of all sizes (up to nearly the system size of 6000 species) are 
present. The results shown are for a system with 6 levels 
and 1000 niches per level, a speciation rate of n — 0.02, and 
a extinction probability (due to environmental changes) of 
p = 0.01. The results are only very weakly dependent on 
the values of the parameters, b Time sequence of extinc- 
tion and origination events. The origination curve is shifted 
downward by 1000 for clarity. Note the strong correlation 
between the two curves, in agreement with empirical obser- 
vations • c Probability density function of events size. 
The distribution is well described by a power law with an ex- 
ponent r = 1.97 ± 0.05, which is consistent with empirical 
measurements 1231. 
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FIG. 3. Fractality of species diversity, a Number of 
species in the model as a function of time. The thicker dot- 
ted line shows the number of species at intervals of 128 time 
steps. The continuous line, shown for a shorter period of time, 
is sampled every time step. Note the complex structure of the 
curve at very small time scales, which suggest that the fluc- 
tuations have a self-afnne |^lj structure |j. b We investigate 
the power spectrum of the signal in a and find that it scales 
as a power law with an exponent /3 = 1.95 ± 0.05, confirming 
the fractal nature of the fluctuations in the number of species 
for the model. 




FIG. 4. Correlations in the fossil record and in the model. 
a We compare the results for the model with the empirical 
data found in |7J. For the model, we consider two sequences, 
one with 512 points (black circles) and another with 4096 
points (dashed line). The figure shows that the scaling be- 
havior found for data and model is similar. We find that for 
about one order of magnitude the data for the shorter se- 
quences appears to scale as a power law with an exponent 
— 1. However, it seems that such scaling does not hold for 
longer sequences, for which the power spectrum becomes flat, 
suggesting that the sequence crosses over to uncorrelated be- 
havior (white noise), b We use detrended fluctuation analysis 
jn] to test the results of the power spectrum. We find F(t), 
which measures fluctuations at different time scales, to scale 
as a power law with an exponent close to 1 for about one 
order of magnitude. In the inset, we show the values of the 
exponent for a local fit to a power law. Again all curves seem 
to behave in similar fashion. However, the results suggest 
that no true scaling regime exist for time scales shorter than 
300; then, the exponent becomes 1/2 which would suggest an 
uncorrelated process. 
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